xy'=y sqrt(x^2-y^2)